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Abstract: Recently a new Lagrangian framework was introduced to describe interactions 
between scalar fields and relativistic perfect fluids. This allows two consistent generaliza¬ 
tions of coupled quintessence models: non-vanishing pressures and a new type of derivative 
interaction. Here the implications of these to the formation of cosmological large-scale 
structure are uncovered at the linear order. The full perturbation equations in the two 
cases are derived in a unified formalism and their Newtonian, quasi-static limit is studied 
analytically. Requiring the absence of an effective sound speed for the coupled dark matter 
fluid restricts the Lagrangian to be a linear function of the matter number density. This 
still leaves new potentially viable classes of both algebraically and derivatively interacting 
models wherein the coupling may impact the background expansion dynamics and imprint 
signatures into the large-scale structure. 
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1 Introduction 

It is generally believed that in order to describe the complete universe history, from its very 
early stages, namely the inflationary epoch or alternatively the bouncing regime, up to the 
late-time accelerated era, one needs to introduce additional degrees of freedom, beyond 
either general relativity or the standard model of particle physics. If these new fields are of 
gravitational origin then we may refer to “modified gravity” [1, 2], which exhibits general 
relativity as a particular limit; and if the extra degrees of freedom belong to the matter 
content of the universe, we have “dark energy” and “dark matter” [3, 4], although one 
can construct combined scenarios of the above extensions, for instance introducing various 
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non-minimal couplings, and typically shift between them, partially or completely, using 
suitable transformations. 

There are numerous proposed realizations of the dark sector, but most of them allow 
at least an effective description in terms of scalar fields. Dark matter is usually modeled 
as a cold (i.e. pressureless) perfect fluid component, and in its simplest and prototypical 
form, dynamical dark energy can be modelled as a “quintessence” [5-7] scalar field rolling 
down its potential, as the inffaton field in the early universe at vastly higher energy scales. 
Generalizations include non-canonical kinetic terms [8-11] and non-minimal couplings to 
e.g. the Ricci curvature [12-15], the Gauss-Bonnet invariant [16-18], torsion and nonmetric- 
ity [19-21], or the derivatives of the field [22-24]. The most general viable action would 
encompass the Horndeski-type [25, 26] and even more general [27, 28] theories. From a 
particle physics point of view, we generically expect the new scalar degree(s) of freedom 
to interact with each other and with (dark) matter [29-42]. Such an interaction could in 
fact potentially alleviate the coincidence problem, namely the comparability of the dark 
matter and dark energy densities today may be more naturally explained if they don’t scale 
independently throughout the evolution of the universe. 

As the microscopic theory of the dark sector is completely undisclosed, it is useful 
to seek a general parameterization of the possible cosmological implications of theoret¬ 
ically consistent models. Recently, a new approach to constructing interacting theories 
of the dark sector was presented in [43, 44]. There, dark matter is described as a per¬ 
fect fluid, using Brown’s Lagrangian formulation of relativistic fluids that employs a set 
of Lagrange multipliers [45]; see Refs. [46, 47] for applications of the same formalism to 
non-minimally coupled theories. By adding a scalar field into the theory and allowing suit¬ 
able non-minimal couplings with the fluid variables, in particular either the fluid’s density 
functional or its four-velocity, one then obtains classes of theories in which the interacting 
dark matter component has unusual properties compared to the previously studied coupled 
quintessence models. We call these new theories the Scalar-Fluid theories. An immediate 
benefit of the Lagrangian formalism is that the system of conservation equations is auto¬ 
matically satisfied, and given the form of the Lagrangian, the interaction terms are uniquely 
determined up to all orders. This is not the case in the approach adopted in e.g. Refs. [48- 
57], where the phenomenological parameterisation of the coupling terms directly in the 
conservation equations has to be done “by hand” at both the levels of background and of 
perturbations^. Nevertheless systematical, general parameterisations of coupled dark sec¬ 
tor cosmologies have been very recently constructed within the frameworks of effective field 
theory [60], pull-back-formalism for fluids [61], and the parameterised post-Friedmannian 
formalism [62]. 

In this paper we investigate the behavior of the scalar-fluid cosmological scenarios at 
the perturbation level. In particular, we uncover the implications of the new couplings 

^Naive phenomenological parameterisations may fail to capture the physics of more fundamental ap¬ 
proaches (as one may argue on general grounds [58] and see directly by comparing with the results derived 
from explicit high energy physics theories [59]); furthermore, unphysical instabilities [51, 52] or even in¬ 
consistencies with the covariant stress energy conservation [32, 48] are possible unwanted artefacts of such 
parameterisations. 
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these models feature, to the large scale structure formation, with the aim of confronting 
the theories with the high-precision data on the galaxy distributions and gravitating matter 
sources in the universe. To begin, we will review the theoretical framework of Refs. [43, 44] 
in section 2. We shall then first present the cosmological equations in Section 3 and in 
Section 4 derive their relevant quasi-static limit for the models with algebraic coupling 
between the dark components, and separately, analyze the qualitatively different scenarios 
with derivative couplings. In both cases we are able to understand the generic features of 
linear perturbation evolution in these classes of models already by our analytic considera¬ 
tions. Conclusions and perspectives that emerge from our analysis are provided in section 
5. 

Some additional material is confined to the three appendices: in A we present the per¬ 
turbation equations in an alternative formulation, in B we sketch a convenient approach 
to solve the perturbations equations numerically making some first steps towards full nu¬ 
merical analysis of the new models, and finally in C we consider, from the point of view of 
high-energy physics motivated couplings, some natural generalisations of the scalar-fluid 
theories. 

2 Scalar-fluid theories 

In this section we briefly review scalar-fluid theories, in which the fluid and the scalar field 
describe the dark matter and dark energy sectors respectively, based on [43, 44]. In the 
first subsection we briefly review the basic formalism, then in 2.1.2 we generalize it to the 
cases of algebraic and derivative couplings, and finally we extract the gravitational field 
equations in 2.2. 

2.1 Lagrangian formulation of a scalar field interacting with a relativistic fluid 

We will employ the Lagrangian formalism for relativistic fluids as presented in Ref. [45], 
since it is complete in the sense of incorporating thermodynamics too. Alternative La¬ 
grangian formulations have been developed by many authors, see for example Ref. [63] for 
an earlier and somewhat simpler formulation and Refs. [64-66] for a more general frame¬ 
work of an effective theory of fluids. 

2.1.1 Minimal theories 

In order to construct a Lagrangian of a scalar field interacting with a relativistic fluid, let 
us first start from the Lagrangian of a minimally coupled relativistic fluid following [45]. 
Such a Lagrangian can be expressed as 

^{minimal) ^ g) + JM + gQ , (2.1) 

where is the metric, with g its determinant. In the above expression p[n, s) is the 
energy density of the matter fluid, assumed to be a function of the particle number density 
n and the entropy density per particle s. Moreover, we need to introduce the Lagrange 
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multipliers y?, 0 and /3 a) with A = 1,2,3, and where aA are the Lagrangian coordinates of 
the fluid. The vector-density particle-number flux obeys the following relations: 


J^ = V^nU^, 


\J\ = , 


n = 


1^1 


( 2 . 2 ) 


with the fluid 4-velocity satisfying = — 1. In summary, the independent dynamical 

variables of the Lagrangian (2.1) are 5 ^^, J^, s, (p, d, Pa and a^, and the corresponding 
variations will give rise to the field equations. 

Variation with respect to the metric gives the energy-momentum tensor 


T/j^U — P Ufj^Uy + 



{Pfiu T UpUy) , 


(2.3) 


which can be rewritten in the usual perfect-fluid form 

-2 5C^ 




if we identify the pressure p as 


= P9^iu + {p + p)Uf,Uu, 


dp 




(2.4) 

(2.5) 


Additionally, variations with respect to the remaining dynamical variables lead to: 


: 

p + s9^fj, + Pa'a'^ — 0 , 

(2.6) 

s : 


(2.7) 

if : 

II 

0 

(2.8) 

9 : 

C/5 

II 

0 

(2.9) 

Pa- 


(2.10) 


(J^/3a),;. = 0 , 

(2.11) 


with 


P = 


P + P 


n 


dp 
dn ’ 


( 2 . 12 ) 

Relations 


the chemical potential, and where a comma denotes a partial derivatives <9^. 

(2.8) and (2.9) respectively correspond to the particle number conservation and the entropy 
exchange constraint, which is consistent with the fact that p and 9 are Lagrange multipliers. 
They can be re-expressed as 


V^(n[/^) = 0, (2.13) 

V^(snV'^) = 0, (2.14) 

where denotes the covariant derivative. Similarly, according to (2.10), the three La¬ 
grange multipliers Pa constrain the fluid’s 4-velocity to be oriented along the flow lines 
with constant . Note that combining (2.7), (2.10) and (2.11) allows us to extract the 
chemical free energy as 

F = p-Ts = . (2-15) 
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Using the above field equations one can easily show that the fluid energy-momentum tensor 
is conserved, namely = 0. 

At this point, let us for concreteness consider some examples of fluids. The most 

important characteristic is the equation of state w, dehned hy p = wp. A barotropic 

fluid with the constant equation of state wq is modelled simply as p{n) ~ A more 

nontrivial equation of state is that of a “Chaplygin gas”, p = —Alp~°^, with some positive 

1 

constants A and a [67, 68]: this corresponds to p = (A + i+“, for a suitable B. 

The so called Cardassian models lead to similar expansion histories without separating 

dark matter and dark energy [69], and the corresponding fluid description would be given 
1 

by p ~ n (1 -|- An~‘^’^) 1 , where A is again a constant of appropriate dimensions and q and 
ly are dimensionless parameters. 

Finally, imposing equations (2.6)-(2.11) into (2.1) one can obtain the on-shell La- 
grangian of the fluid as [43] 

^M = -'/^P = V^P (on-shell), (2.16) 

where the last equality holds up to total derivatives and only once the equations of motion 
of the non-interacting fluid have been taken into account. Thus, this will not hold in general 
in the following, where we will consider non-minimal couplings of the fluid variables^. 


2.1.2 Non-minimal theories 


Having obtained the Lagrangian for the dark matter fluid, we can now proceed to the con¬ 
struction of the total action of scalar-fluid theories, which includes an interaction between 
the matter fluid and the scalar field that describe the dark energy sector. In particular, 
the total action consists of three pieces: 


5 



+ l^M + 


(2.17) 


In the above expression Cm is the perfect-fluid Lagrangian (2.1) that is now generalized to 
include interactions with the scalar field, £grav is the standard Einstein-Hilbert Lagrangian 


p _ 3 p 

-C-grav — 2^2 ^ 


(2.18) 


with R the Ricci scalar and k the gravitational coupling, related to the Newton’s constant 
G through k? = SttG, and where C^ is the free scalar-held Lagrangian, which we take here 
for simplicity to correspond to a canonical scalar held cj) with a potential V{4>), namely 




-V^g 


]^d^(l)d^(l) + V{4>) 


(2.19) 


In the total action (2.17) we have included the interactions in the matter Lagrangian Cm^ 
which thus incorporates the coupling of the huid and scalar degrees of freedom. The energy 

^In the context of non-minimal curvature-matter couplings [70-72], the above observation clarifies the 
pseudo-issue of the choice of the matter Lagrangian; see discussions and references in section 2.2 of the 
review [72]. 
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density p of matter for example corresponds then to the total energy density including 
the contribution from the interaction with the field. In principle the separation into the 
fluid and field terms is arbitrary in the coupled case [58], but lumping all the matter- 
coupled terms besides the canonical scalar terms in (2.19) to the “fluid sector” appears 
both physically intuitive and technically more convenient^. Moreover, for clarity, we will 
focus separately on the two distinct generalizations (though in principle both could be 
present simultaneously): 

• Algebraic coupling [43]: 

Cm =p{n,s,(l)) + J. ( 2 . 20 ) 


• Derivative coupling [44]: 

Cm = -y/^p{n, s) -b fin, s, + J. (2.21) 

Here J is a shorthand notation for the second term in the RHS of (2.1). In both (2.20) 
and ( 2 . 21 ), f[n,s,<j)) is an arbitrary function of the scalar field </>, of the fluid particle 
number density n and of the entropy density per fluid particle s, and thus one has an 
implicit coupling between the fluid dynamical degrees of freedom (i.e excluding the various 
Lagrange multipliers) and the scalar field. The difference between (2.20) and (2.21) is that 
the former represents a general coupling in the absence of derivatives, while the latter is 
a coupling that one can construct in the presence of scalar-field derivatives, which should 
then need to be coupled with the the fluid vector-density particle number flux J^. 

In summary, the total action (2.17) describes theories where the new scalar-fluid in¬ 
teractions stem from the new field-dependent terms ( 2 . 20 ) and ( 2 . 21 ) implemented in the 
Lagrangian variational approach to relativistic fluids. The functions p and / in the ac¬ 
tion determine the coupling completely, the system is automatically consistent with the 
Bianchi identities, and there are no ambiguities to be fixed by hand, unlike in models 
based on phenomenological parameterizations of the conservation equations. 

2.2 Field equations 

Let us now use the total action (2.17) of scalar-fluid theories, in order to extract the held 
equations, by varying with respect to , cj), J^, s, p, 9, f^A and We will do this 
in a unihed notation and specify the relevant terms in the equations of motion separately 
for algebraic couplings [43], i.e for the interacting Lagrangians (2.20), and for derivative 
conplings [44], i.e for the interacting Lagrangians (2.21). 

In the analysis below, we will neglect the role of the entropy density s. In fact from 
equations (2.13) and (2.14), one hnds that in general 

9^5 = 0 , ( 2 . 22 ) 

®We present the equations in an alternative formalism that isolates the coupling terms as Cm = 
^{minimal) ^ appendix A. This coincides with the notation originally used in [43, 44]. In the 

following we will however use the simpler notation of Eq. (2.17). 
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a relation that holds regardless of the couplings. This implies that the entropy density is 
always a constant and we have an adiabatic process. Did one want to include entropic 
effects, one needed to modify the constraints by adding dynamics to s in the Lagrangian. 
However, as we focus here on the effects of the couplings we will for simplicity neglect the 
possible intrinsic entropy of the fluid. Hence, the functions coupling the scalar to the fluid 
in (2.20) and (2.21) can be considered to depend only on n and (p in what follows. 

The field equations given by the variation with respect to the metric, retain their usual 
form, namely 

+ rW) , (2.23) 

where we recall from (2.4) that 


_ -2 SCm f , \ 

= 7^^ = (P + P) > 

and similarly for the scalar field, for which the canonical Lagrangian (2.19) gives 


(2.24) 


^ 


du4> - g^u 




(2.25) 


The crucial issue is definitely that the components of the energy momentum tensor 
(where for simplicity we have omitted the superscript M for the matter fluid), depend 
also upon the scalar field (and possibly its first derivatives). Consequently, the individual 
energy momentum tensors are not conserved, though their sum, due to the diffeomorphism 
invariance of the matter action and in accordance with the Bianchi identity = 0, 

is, namely 

v^(r^, + rW) =0. (2.26) 

Following the usual notation, we can parametrize the interaction with a four-vector as 


= Qy = -V^^T^y . 


(2.27) 


It remains to deduce the non-conservation terms in the two cases we are considering. 
We obtain, respectively, the following results: 


Algebraic coupling (Eq. (2.20)): 


P = p{nA), 
dp{n, p) 


p = n- 


dn 


- P{n,4 >), 


^ _ dp{n, (p) ^ ^ 
Qu - o i * 

ocj) 


Derivative coupling (Eq. (2.21)): 

P = pin), 
dp{n) 


p = n- 


Qy = -n 


dn 
2d f{n,p) 
dn 


( ^ 2dfin,p)^^xr, j, 
-p{n)-n - —U dxp, 


dn 


^xU^dyP. 


(2.28) 


(2.29) 
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In both cases, the Klein-Gordon equation follows from (2.27) and (2.25) as 


(□</. - V') d^cj) = . (2.30) 

The stress energy conservation equations follow from (2.27) and (2.4), and can be written 
as 


+ {p + p) , (2.31) 

(p + p) + V^p = - (p^^ + , (2.32) 

when projected along the fluid flow and orthogonal to it, respectively. Our notation slightly 
differs from Refs. [43, 44], however the results coincide, and we refer the reader to those 
works for more details on the derivations. 

3 Cosmological equations 

In the previous section we presented the formulation of scalar-fluid theories, in which one 
incorporates the scalar-fluid interaction, i.e. the interaction between dark energy and dark 
matter, straightaway at the Lagrangian level, instead of inserting it by hand at the level 
of equations of motion. We extracted the covariant field equations, which can be used in 
the relevant applications. Here we specialize to cosmology, first deriving the equations of 
motion adapted to the background describing the overall cosmic expansion, and then to 
linear perturbations, describing the formation of cosmological large-scale structure. 

3.1 Background 

In particular, in order to explore the cosmological applications of the above new class of 
theories, we focus on the Friedmann-Robertson-Walker (FRW) metric of the form 

ds^ =-dt^-h a(t)^ -k , (3.1) 

with a{t) the scale factor, K = —1,0,1 for spatially open, flat or closed geometry respec¬ 
tively, and with dfl^ the two-dimensional line element of a sphere. Additionally, all the 
variables that appear in the equations of the previous sections are assumed to be homoge¬ 
neous, that is depending only on the cosmic time t. Finally, concerning the fluid 4-velocity, 
in comoving coordinates it becomes = (1,0,0,0). 

In the FRW geometry, the Friedmann equations of the scalar-fluid theories read [43, 44] 

3— -|- ^p -|- ^ 

- , (3.3) 

where the dot denotes derivative with respect to t and we have introduced the Hubble 
parameter H = a/a. As expected, these equations retain their usual form. Moreover, we 
obtain the conservation equation for the number density as 


^ + 2H + 3H^ = 


P+7;4>‘^ 


h + 3Hn = 0 ^ n ^ a 


(3.4) 



Had we kept the entropy s in the system, we would have recovered its conservation equation 
at the FRW background as just s = 0. Finally, the conservation equations feature the 
couplings as follows: 

p + 3H{p + p) = Qo, (3.5) 

^ + 3H<j) + V'= -Qo/(j). (3.6) 

In these conservation equations the scalar-fluid interaction is taken into account by the 
time-component of the coupling vector Qfj^, namely Qo- In th® two cases we consider, the 
coupling turns out as: 

• Algebraic coupling (Eqs. (2.20),(2.28)): 

Q« = p. (3.7) 

• Derivative coupling (Eqs. (2.21),(2.29)): 

Qo = -3Hn^^^. (3.8) 

an 

The cosmological background dynamics of scalar-fluid theories at the background level were 
studied in terms of phase-space analysis in [43], for the algebraic coupling case, and in [44] 
for the derivative coupling case. The results in those two works imply the existence of 
many new realizations of viable cosmological evolutions in the scalar-fluid theories. Hence, 
it would be interesting to understand their implications to the structure formation. 


3.2 Scalar perturbations 

In the present work our main focus is on the behavior of the scalar-fluid cosmological 
scenarios at the perturbation level. We will first extract the general equations for scalar 
perturbation, and then use them in order to examine the signatures of these models in the 
structure formation. We shall continue in a unified notation until the end of this section 
and then perform a detailed analysis separately for the algebraic couplings and for the 
derivative couplings. 

We will work in the Newtonian gauge, unless otherwise specified. The Newtonian 
frame is also called the longitudinal gauge, since there the perturbed metric can be written 
(in Cartesian coordinates) as 

ds2 = -(1 + 2$)dt2 +- (l-2T)a^(t) -^ ^ ^ 2 ^ ^^2^ (3 9) 

[l + iR(x2+y2 + ^2)]2 I I’ 


where T, <I> are functions of all coordinates. Since in our case all the matter sources can 
be considered as perfect fluids, no anisotropic stresses appear in the considered scalar-fluid 
models. Hence, the off-diagonal ij-components of the Einstein field equations always read 
as 


didj 


l + ^^K{x^ + y^ + z^) 


(T-4>) I = 0 


(3.10) 
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and hence we can immediately deduce that 


$ = ^', (3.11) 

as expected in the absence of source imperfections (i.e. off-diagonal pressures). Thus, in 
what follows we will simplify the equations by setting equal to T. 

Let then first determine the matter variables that are going to be perturbed. In scalar- 
fluid theories the fundamental matter variables of the coupled system are^ n, (j) and {s 
is neglected due to adiabaticity as we mentioned above) and therefore we should consider 
the perturbations of these quantities as 


4>t-^4> + d4>, n n + 6n, Ufj_ Ufj_ + 6Ufj,, (3-12) 


where (j), n and U^j, are now the background quantities, and 


5U^ = {--^,div) , 


(3.13) 


with V the scalar perturbation of the matter fluid velocity®. We parametrize the fluctuations 
in the energy density as d = 5p/p, and for the pressure perturbation we simply write 6p. 
These two quantities have a different relation to 6n depending on the perturbation of the 
general definitions of p and p given in Eqs. (2.28) and (2.29). Below we will also give the 
relation between 5p and 6p in detail for the two classes of models studied here. 

We are then ready to write down the remaining field equations. There are three 
independent Einstein equations remaining, since we have already used one to set $ = T. 
These can be written as follows, lumping the matter and field fluctuations on the RHS and 
the geometric fluctuations on the LHS: 


K e 


6^ - ^ - kW T - 3H^ 


+ ^5(j) + V'6^ 




2 1 


K 1 


T 4iLT + {2H + 3H^ - ^ 1 ^ = 


K 

T 


{p + p)v - 
[Sp + 4>5cj) — 


(3.14) 

(3.15) 

(3.16) 


Here k is the wavenumber of the fluctuation (where in A:-space with suitable 

scalar harmonic eigenfunctions that for K = 0 reduce to plane waves). To complete the 
system, we need also the conservation equations for the fluid and the field. We will study 
these separately in the two cases at hand. 


4 Structure formation 

In this section we explicitly elaborate the way from the general perturbation equations to 
the physics of structure formation at the linear regime. 

"^The actual variable upon which the scalar-fluid action is varied, is and not , however according 
to (2.2) these are related through = y/—gnU'^, and thus one can always consider instead of J''. 
®The fact that 5Uo = — comes from perturbing the relation = — 1. 
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4.1 Algebraic couplings 

In this subsection we analyze the scalar perturbations in the case of algebraic couplings 
[43], characterized by Eqs. (2.20), (2.28) and (3.7). We will need various derivatives of the 
coupling functions, and it turns out to be useful to define the following auxiliary background 
quantities (all with the dimensions of 1/M): 


X = 


{l + w) p' 


y 


_ p,<t> 

~ p 




p 


U=' 


(l + w) p' 


r = a . (4.1) 


(1 + w)p ' 


Furthermore, we introduce the quantity that is the sound speed square of the fluid in the 
rest frame of the field: 


„2 _ „ P,r. 


p,r 


(4.2) 


If the field is decoupled from the fluid, this coincides with the usual definition of the 
adiabatic sound speed squared (the background quantity) = p/p, and at small scales 
it coincides with the usual definition of the sound speed squared (perturbation quantity) 
that is evaluated in the rest frame of the fluid, = dpfdp, where the hat indicates that the 
fluctuations are considered in the comoving matter gauge wherein D = 0. It turns out that 
in the limit we will consider, the effective sound speed is indeed given by Eq. (4.2), since 
at sub-horizon scales the difference between the two rest frames is negligible. However, the 
reader is warned that fundamentally our definition is slightly different from the conventional 
one, for purposes of convenience in the present case. At the background the evolution of 
the pressure is given by (recall that the equation of state is defined by re = p/p) 


ri; = (1 -|- tc) 3i7 {w — cf) + {x — y)f 
and for the perturbations the total fluid pressure is given as 

— = cl5+ [(1 -h u;) X - (l cl) y] df , 


(4.3) 


(4.4) 


which shows that (4.2) would indeed give the total pressure perturbation of the fluid only 
in the rest frame of the field df = 0, and that in general (4.2) differs from c^. 

We can then write the conservation equations for our coupled scalar-fluid system. The 
Klein-Gordon equation for the fluctuation of the field is 
ri.2 -I 

d<p + 3Hdf)+ ^ + V" + p{z‘^ -xy) d(f + 2{py + V')'i’-Af)'^ + pxd = 0. (4.5) 

The continuity and the Euler equation for the fluctuations in the matter fluid come out 
respectively as 


<5 + 


and 


3H [cl — w) + (j) {y — x) (5 — (1 -k rc) ( —v -k 34' j = yd(j) 




V — I 3Hcl — 


+ ^3H [(l + cl) y - {l + w)x] -<j) [xy - z^) j d(p , (4.6) 

(4.7) 


u -k T = 


1 -k rc 


y- X 


1 -k rc 
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This completes the perturbation system for the models defined by Eq. (2.20). Definitely, 
when we switch off the coupling by setting x,y,z 0, we recover the usual conservation 
equations. In the extensively studied case of conformally coupled cold dark matter, cor¬ 
responding to yO = e^^niQU, we have simply x = y = z = u = r = l3, and the equations 
simplify to 


C-coupled CDM : 


' 5 ^ + 3H5^ + (§ + V”) + 2E'T - 40T 

< 6-^v-34> = P6<P, 

— f3^v -|- T = —/36(j) ■ 


-f3p{6 + 2^) , 

(4.8) 


Obviously, allowing /? to be a general function of both the scalar field and the number 
density of matter particles results in a richer structure of the equations. We proceed now 
to explore their cosmological implications. 


4.1.1 The quasi-static limit 

Cosmological large scale structure can be probed most efficiently at sub-horizon scales, since 
beyond that we face problems related to both systematic difficulties and cosmic variance. 
At the subhorizon scales, we can considerably simplify the system of perturbation equations 
by considering taking advantage of the so called quasi-static limit, that will be explained 
in the process of exploiting it below. 

To begin, we combine equations (3.14) and (3.15) to obtain the generalized Poisson 
equation in the quasi-static limit: 


6 ^- 


+ 


(p + P + <))2) 


T = 


p6 (v' + 54> + 4>S<j) 





(4.9) 


In the first line we have eliminated the T and at the same time obtained the density 
perturbation in the comoving matter gauge, 5 = 6 — 3H{1 + w)v, which is the quantity 
we are interested in®, as we presumably are comoving with matter, and thus this is the 
gauge corresponding to physical observables such as the matter power spectrum, ~ 

Let us now discuss in detail how we can justify, within our approximation, to drop each 
of the additional terms that appear in the hrst line of the equation above. Firstly, at the 
sub-horizon scales we consider, k‘^ ^ K, as the spatial curvature scale of the universe is 

®There is also a subtlety involved here. In fact we could argue that the quasi-static limit of (3.14) 
would directly yield the Poisson equation (4.9) with S replaced by the Newtonian gauge S. However, in 
the presence of non-negligible pressures, the velocity fluctuations and thus the difference in the gauges can 
be significant, and though not naively apparent from (3.14) alone, the quasi-static limit of the full system 
would be inconsistent if we removed the hat from (4.9). Care has to be taken with the derivatives as it 
turns out that in fact 4' ~ ti^c^pS/2 in the subhorizon approximation (and in this sense the quasi-static 
approximation is not strictly valid), though still 4* ~ 0{H’i!) and thus we are able to close the system 
consistently. In particular, the approximated system of equations we eventually obtain, described by the 
result (4.11), satisfies all the field equations and the Bianchi identities. 
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constrained by observations to be at least beyond the present horizon. Secondly, since 

^ well as K^V,K^(jP‘ < and, by definition of the subhorizon limit, we 

consider scales for which k'^ S> the gradient term in the LHS square brackets clearly 

dominates over the other terms. Similarly, in the RHS the matter perturbation (term ~ 6) 
will dominate over the held fluctuation (terms ~ 64>,5(j)). At small scales the held may 
perform rapidly small oscillations, and thus we only need to take into account its averaged 
amplitude. This we will soon see to be given by the coupling to matter, in such a way 
that the held huctuations are indeed suppressed in comparison to the matter over densities, 
k6^ ~ {H/kf5. 

Therefore, using (4.9), we can relate the gravitational potential 'k directly to the matter 
overdensity 6, and turn the evolution equation for 4', the space-space component (3.16) of 
the held equations, into an evolution equation for 6. However, in order to close the system 
we also need to express the scalar held perturbation (and its derivative) in terms of 6. This 
can be accomplished by using the quasi-static limit of the Klein-Gordon equation (4.5), 
namely 

5(j) = xp5. (4-10) 

We arrived at this simple relation by repeating the arguments above that we used to justify 
Eq. (4.9). The time-derivative of the gravitational potential is of the order of magnitude 
'k ~ ff'k, since, as in general also for other huctuations, the time-scale of the evolution is 
given by the cosmological expansion rate: d/dt ~ H. We have then everything needed to 
close the system of equations. In equation (3.16), where the pressure perturbation is now 
given by (4.4), we hrst substitute 4* and its derivatives using (4.9) and 5(j) and its derivatives 
using (4.10), and then we elaborate the resulting equation using the background formulas 
(3.2)-(3.6) and (4.1)-(4.3), obtaining 



5 + 


(2 — 6 u; -|- 3c^) H -|- 2y^ 


s + ci 




+ CiH<t) + ) <5 


(4.11) 


where 


Co = 


K 


(1 — 2>w) {l + w)p+ [4cg — 3w) -k (l5rc — 90^) + {2 + 3rc) ^ ^ 

-k [(1 + w)x^ -{l + cfj xy + y'^]p-{x- y) V ', (4.12) 

Cl = (5 -k 6 c^ -k 3rc) x — (5 -k 3c^) y — 3r , (4-13) 

(4.14) 


C 2 = V? — x^ — xy — Ay^ . 


This result is applied in Appendix B.l to analyze a specihc class of Lagrangians. 

We close this subsection by checking some special limits of this equation. In the 
following let us set K = 0 for simplicity. Firstly, if we switch off the scalar held, we recover 
the exact evolution equation for overdensities in the uncoupled matter distribution: 


Matter only: 


5 + (2 - 6u; + 3 c 2) i7<5 + ^ [1 + (8 - 3u;) u; - 6 cf] H^5 



(4.15) 


- 13 - 






The hatted sound speed should be used in the entropic case (though not considered here). 
The mass-varying cold dark matter, described by equations (4.8), results in 

C-coupled CDM: 6 + ( 2 H — S = p6. (4.16) 

In addition to the friction term on the LHS, on the RHS the coupling always adds to the 
effective Newton’s constant: the fifth force between CDM particles that is mediated by the 
scalar field is always attractive^. Both of these effects are signihcant in the perturbation 
evolution, and generically modify the matter power spectrum too much to comply with 
observations, in cases when the coupling (3 is large enough to qualitatively change the 
background evolution: this is a basic difficulty in addressing the coincidence problem within 
the standard coupled quintessence models featuring a canonical scalar field coupled to cold 
dark matter [32]. 

In general, however, modifications arising in the form of the effective sound speed term 
in (4.11) are still more disastrous to the matter power spectrum due to their fe^-enhanced 
influence at smaller scales. Therefore, we can exclude any appreciable nonlinear dependence 
of the fluid Lagrangian upon the number density. We will return to discuss this result in 
Section 5, after noting that it applies also for the derivatively coupled models in section 
4.2. 

4.2 Derivative couplings 


Let us now consider the case of derivative couplings [44], characterized by Eqs. (2.21), 
(2.29) and (3.8). In order to simplify the notation in what follows we define 


X = n 


dcpdn ’ 


2 5/ 

y = ^ ^ 

an 


z = n 


d’n? 


(4.17) 


From (2.21), we now obtain for the pressure at the background and at the perturbation 
levels 


P = np^n - p-y^. 


ip 


^ \ 5 + y(t>^ — x(j)S4> — y5(j >, 


I+ Wq 

respectively, where 5 = 5p/p, wq (the “bare” EoS parameter) is given by 


(4.18) 

(4.19) 


_ P,n ^ y<P / . 

Wo — n -1 = tc -I-, (4.20) 

P P 

and Cg is defined exactly as in Eq. (4.2), although now p does not depend on the scalar 
field. In analogy with wq, corresponds now to the “bare” sound speed and we remind 
that it is not the “full” as explained after (4.2). The evolution of the equation of state 
is given by 


w = 2>H ( 


(1 + rco) (re - c^) - ^ + 377 {2y + z) (p - y"^ EV} • (4-21) 


^This result holds with generalized “disformal” couplings too [73], though could be avoided at least with 
a (pathological) phantom scalar field [74]. 
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It is a curious property of the coupling that it adds only a contribution to the effective 
pressure. This allows for example to construct a de Sitter expansion dominated by the 
canonical kinetic energy of the scalar field: the scalar can have a constant time-derivative 
due to the coupling [44] . 

The perturbed scalar field equation reads 

S(j) + 3H6(j) + - ^Hx -t Scj) - {3Hy - 2V') 4- - - Sy) = 3H ^ ~ ’ 


(4.22) 

In contrast to the algebraically conpled models, here the scalar flnctuations are sourced 
also by the velocity perturbation in the matter field. The continuity equation is simply 


S + 3(Cg 


wo)H6 - (1 -h Wo) 



v + 3'I> 


= 0 , 


(4.23) 


and thus retains its nsual form unaffected by the coupling. This can be seen from the 
covariant equation (2.31), that for the derivative couplings (3.8) reduces to (1 -|- 

wo)pV = 0 , where the scalar-field dependence has vanished identically because the 
“fifth force” is orthogonal to the fluid flow as shown in [44]. In the Euler equation, which 
is (2.32) adapted to cosmology, we however obtain nontrivial conpling terms as the spatial 
components are subject to the interaction: 


{l+w) pi] + I 




3H 


cIp 


(1 -h Wo) cIp -{2y + z)^-y^ + - yE'j v 

- (5 -I- (1 -I- Wo) = {^Hy + 6(t> + yS^. 


(4.24) 


Note that in the equations z appears only in the combination 2y -|- z. 


4.2.1 The quasi-static limit 

We can now consider the derivatively coupled models in the quasi-static approximation. 
Since onr focns is upon the novel effects arising from the gradient-type interaction, we 
consider the matter source to be cold, p oc n by itself for simplicity, that is we set wo = 
c2=0. 

In addition, we restrict the analysis to the special class of couplings for which / oc 1/n. 
The reason is that in the first place the quasi-static approximation would not be justified 
in the presence of any nonlinear function n^f; such a term would contribute an effective 
pressure as a source in complete analogy with a nonlinear p[n) term, the difference now 
being that the pressure term due to the coupling / appears with the gradient-type prefactor 
and it is thus strongly enhanced at small scales. The two effects combined, extra pressure 
and extra gradient, would probably lead to disaster for structure formation at subhorizon 
scales (nnless the coupling is tuned to be extremely tiny), and since then the scalar field 
would cluster together with matter also inside the horizon, the quasi-static approximation 
would not be valid for these models®. 

®In fact, if we follow the steps of subsection (4.1.1) in such a case, we would find that the solution for the 
50-derivatives, in addition to derivatives of 4, are of the order of 5, and then we could not close the system 
consistently (one could though reduce it to two coupled second-order equations without the quasi-static 
approximation) - see also the comments in the previous case in the footnote 6 of section 4.1.1. 
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Hence, in what follows we specialize to cold dark matter with the derivative coupling 
z = —y/2 and x = y', where a prime denotes differentiation with respect to 4> and y{4>) has 
to be taken as a general function of the scalar field. This simplifies the system considerably. 
Similarly to the previous case, the time-time component (3.14) of the field equations reduces 
to the usual Poisson equation in the quasi-static limit: 



(4.25) 


Thus, we can again relate the gravitational potential ^ directly to the matter overdensity 
6, and turn the evolution equation for the space-space component (3.16) of the field 
equations, into an evolution equation for 6. However, to close the system we also need to 
express the field perturbation (and its derivative) in terms of 6 (and its derivatives): 


64 ) 


= -yv = -y 



(4.26) 


where we used the appropriate limit of Eq. (4.22) for the first, and of Eq. (4.23) for the 
second equality. Then, after some algebra we arrive at the following closed form equation 
for (5: 


5 + 


2H + 


y'y + y' {‘^y- 


p + y[y- 


i 


1 + 


y\y 


~i 


(4.27) 


We stress that now the effective gravitational constant can receive positive as well as 
negative contribution from the coupling, whereas in the most standard cases the fifth force 
is always attractive, as is clear from (4.16). 

In appendix B.2 we briefly apply the above results to a concrete example, considering 
the models specified by / = —^77o/(^^) (so for which y = 4.Ho/k is just a constant), whose 
background phase-space analysis was performed in [44]. 


5 Discussion and outlook 

In this work we studied the formation of cosmological large-scale structure in the new 
class of Scalar-Eluid theories. These theories are based on an extension of the Lagrangian 
formalism of effective perfect fluids, that allows for new scalar-field-dependent terms in the 
fluid Lagrangian, introducing non-minimal interactions of the perfect fluid and the scalar 
field. We explored the cosmological implications of the new couplings of both the algebraic 
type (i.e. when the fluid energy density may depend non-trivially both upon the number 
density of particles and the scalar field, thus generalizing the usual Yukawa-type coupling 
p = m{(p)n —>• p{(p,n)) and of the derivative type (i.e. where the interaction is between 
the four-gradient of the field and the four-velocity of the fluid, contributing only to the 
effective pressure term, p —)• p — 

The cosmological equations were derived up to linear order in perturbations, and fur¬ 
thermore the perturbation system was reduced to a simple second order closed evolution 
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(master) equation for the observable matter structures S at the relevant subhorizon limit, 
where we were able to employ (partially) the Newtonian quasi-static approximation. The 
main results of the investigation were given as Eq. (4.11) for the algebraic case and as 
Eq. (4.27) for the derivative case. Compared to the standard result of Eq. (4.15), we 
observe three different modifications in the matter perturbation evolution: 

• Additional friction terms. In addition to the usual Hubble-friction that slows down 
the clumping of matter due to the expansion of the universe, the couplings typically 
introduce additional friction terms that are given by the evolution of the scalar field. 

• Modifications of the effective Newton’s constant. Matter falls into the potential wells 
induced by matter, but the rate at which this occurs can be modified by the coupling. 
This can be interpreted as an extra force between matter particles mediated by the 
scalar field, and described as a varying amplitude of the effective Newton’s constant. 

• Effective pressure perturbations. Any nonlinear dependence on the number density 
of the fluid particles in the Lagrangian means nontrivial self-interactions. These 
will introduce pressures in the fluid, thus implying propagation of sound waves. In 
the matter fluctuation equations these appear as gradient-like terms that become 
increasingly important at small scales. 

Given any specific model, i.e. given any specific Scalar-Fluid Lagrangian, these effects can 
be straightforwardly determined from equation (4.11) or equation (4.27), which can then 
be conveniently integrated to obtain the growth rate or the matter power spectrum as a 
function of the redshift; see appendix B. 

In the context of interacting dark matter - dark energy cosmologies, the sound speed 
term is the most dangerous for the viability of the models. Such terms have been known 
to appear in coupled three-form models [24] and in unified dark matter scenarios (such as 
the Chaplygin gas or Cardassian expansion models [67-69]), and it is the culprit which 
makes these constructions ruled out - or rather constrained to a tiny parameter region 
where they are completely indistinguishable from the standard ACDM model, at least 
by any realistic background measurement [75, 76]. For some modified gravity models 
that generate such effective pressure perturbations as source for matter fluctuations, the 
perturbation constraints have been also shown to tighten the viability bounds ensuing from 
background expansion alone by several orders of magnitude [77, 78]. We therefore know 
that any nonlinearities in matter density, either in the minimally or in the coupled matter 
terms in the Lagrangian, could not have any appreciable influence on the background 
evolution®. 

Even restricting to models which exhibit no effective pressure gradients at the per¬ 
turbation level, we have at hand new classes of potentially viable interacting dark sector 
cosmologies. The constraints we can impose on them by a detailed comparison of the the- 

possible caveat is that including both non-standard properties for dark matter and non-trivial in¬ 
teractions, the new effects could in principle cancel each other in order to eliminate the dangerous sound 
speed term [79] . This would however appear to require a very fine-tuned conspiracy. 
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ory predictions and the cosmological data remain to be studied. The present work provides 
the tools to carry out such analysis up to linear order in cosmological perturbations. 

Let us end with the remark that there might be interesting scalar-fluid theories besides 
the two classes focused upon here. In particular, as discussed in appendix C, high energy 
physics considerations naturally suggest more general forms of derivative couplings than 
the linear one studied here, and in addition one could contemplate also on modifying the 
Lagrangian constraints included in the theory. In this way one could change for example 
explicitly the number conservation of particles, implying scalar-field dependent particle 
creation, which could be phenomenologically useful e.g. for the reheating process after 
inflation. Furthermore, including dynamics for entropic degrees of freedom and extending 
to imperfect fluids, are amongst other possible generalizations to consider in the future. 
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A Alternative representation 


In this appendix we will present the cosmological perturbation equations in the notation 
originally adopted in [43, 44], where the fluid Lagrangian is separated into its minimal and 
interacting part. Action (2.17) is written in this notation as 


5 


/ id I , /I , ^(minimal) 

d X Tgrav +Ccj>+J^M 


+ -Cint 


(A.l) 


where Tgravj and respectively defined by Eqs. (2.18), (2.19) and (2.1), 

while the interacting Lagrangian assumes the forms 


Tint = ~/(n, s, <(>) for algebraic couplings , 
^int = f{n, s, 4>)J>^dfj,4> for derivative couplings . 


(A.2) 

(A.3) 


Note that the meaning of /(n, s, (p) is different for the two types of models and one must 
not be confused by this small abuse of notation. We will refer to action (2.17) as the 
“implicit” representation, while action (A.l) will be called the “explicit” representation. 
Such names derive from the fact that in (2.17) the scalar-fluid coupling implicitly appears 
within the fluid Lagrangian, while in (A.l) it explicitly appears as a separate term. 

The main difference between the two representations is given by the definitions of the 
fluid energy density and pressure. In the implicit representation we have the expressions 
given by Eqs. (2.28) and (2.29), while in the explicit representation the energy density and 
pressure are solely given by the minimal fluid Lagrangian as 


p = p{n, s) and p = n 


dp{n,s) 
dn 


- Pin,s). 


(A.4) 


They thus never contain any dependence on the scalar field, which will then always appear 
in separate terms within the field equations. The relation between the energy density and 
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pressure in the two representations coincides with the “tilde” transformations considered 
in [43, 44] and then generalized to phenomenological models of interacting dark energy in 

[58]. 

In this appendix we will not present the covariant field equations as well as the back¬ 
ground cosmological equations in the explicit representation since these have already been 
given in [43, 44], We will only provide the corresponding perturbation equations following 
the same notation and conventions we adopted in the main text. Again because of the 
absence of anisotropic stresses the constraint 'k = will apply, allowing us to simplify the 
equations. 

For the algebraically coupled models (A.2) the perturbed equations in Newtonian gauge 
are then 







- K^{p + V + f) 



/ df n 
\dn p + p 


+ 1 



6(j) 


5p 

- = 0 , 


(A.5) 


^2 

'if+ H'i> + Y 


p + p + n 


dn 


V - Y^^^ = 0 , 


(A.6) 


^ K^d'^f n? 

2 ^ 2 dv? p + p 


6p + 


fdf 


— n 


d^f 


d(j) dnd(j) 


+ V']6^ - 


+ (2H + 3H^ + —4)^ - 


K 


^ + ^ = 0 , 




a? dcj)^ 


d^f 


n 


6(j) + 3H6(p +{'^ + ^ + V" ]6(j)+ " " 6p -2((f) + 3H(j) ] ^ = 0, 


dnd(j) p + p 

e 


5p + 3H {6p + 6p) - {p + p)^v - 3(p -kp)'!' = 0 


(A.7) 

(A.8) 

(A.9) 


dn dn 


V + 


dndcj) 


— 3Hn 


d'^f d‘^p 
dn^ dn^ 


V + 


d^f 


n 


Sp + 


d^f 


dn^ p + p dnd(j) 


64> 


1 


n 


df , dp 


+r^7^+hir + iw = (A-io) 


dn dn 


where 6p and 5p are related to 6n through the perturbation of Eqs. (A.4), and p, p, 
(j) and n now denote background quantities. Eqs. (A.5)-(A.7) arise from the perturbed 
Einstein field equations, while Eq. (A.8) is derived from the scalar-field equation, and 
Eqs. (A.9) and (A. 10) correspond to the perturbed matter conservation equations. Note 
how the interacting function / explicitly appears in various terms, when these equations 
are compared to Eqs. (3.14)-(3.16) and (4.5)-(4.7). 

Eor the derivatively coupled models (A.3) the perturbed equations in Newtonian gauge 
are instead 

^ - K^p - 'k - 3H'if - ^5p - ’^V'Scf) - = 0, (A.11) 

\ J 2 2 2 
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^2 / ■ df\ K^- 

il + + — {p + p - ) V - —^6(p 


(A.12) 


^ + 4i?^ + 


^ tt2 tt '^^5 2 2 ; ^/ 

-^ + 2i? + —- —r?4>T^ 


2 " 2 
-y«P+y(» 


9n 


^ -J df d‘^f\ n2 


Sp 


+ v']54> + 


K- I 2 9 / 


n-—-4i]6(i) = 0, (A.13) 


^TT ! d^f\ V? 


5?i 


+ ( 2(/) + 6H^ — 3Hn 

e 

-.2 


dn 


'I' + 4(/) - 3n' 


5n 


4' 


- Kn^^v + (-K + 3Hr?Srl- - V" ) Sep - 3HSp - Sp = 0, (A.14) 


dpdn 


p + p 


+ 


\ k‘^ 

3HSp + 3HSpj --^v-34/ 


- 34- = 0 , 


(A.15) 


dp -df 
dn "^^dn 


V + n 


an on^ 


df d'^f 


n 


I T 

-Sp + —^ - n 
p + p dn 


nrr^'^P l^f 
dn^ dn 


+ iSV\ 

on 


12 


d^f 

dpdn 

df 


V H —Sp 
n 


dpdn J 


Sp — n-+^Sp = 0. (A.16) 
dn 


Note that Eqs. (A.11) and (A.15) are not modified by the coupling. This happens because 
all terms arising from the interactions are orthogonal to the fluid flow, as already pointed 
out in [44]. Again the differences between Eqs. (A.11)-(A.16) and Eqs. (3.14)-(3.16) and 
(4.22)-(4.24) are in the explicit terms depending on the coupling, though for the models 
(A.3) the dehnitions of p in the implicit and explicit representations coincide. 


B Specific solutions 

In this appendix we make some first steps towards a numerical elaboration of structure 
formation in the coupled models. We sketch an approach to implement background phase- 
space variables in order to conveniently solve the main equations (4.11) and (4.27) for the 
algebraically and derivatively coupled models respectively, and on the way we obtain some 
explicit solutions in a couple of special cases where the background admits an analytic 
(fixed point) solution. The systematic analysis of the full cosmological solutions in realistic 
models and their confrontation with the data is left for a future project. 

B.l Algebraic couplings 

Consider the model specified by 

pin, P) = An^+^° + ^ . (B.l) 

This is equivalent to the model proposed in [43], where in addition to the “bare” energy 
density described by the equation of state wq , we have a coupling term proportional to the 
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power a of the “bare” energy density and the exponent of the field with slope /3. We can 
parametrize the background evolution in terms of the phase-space variables^^ a, X, Y and 
Z, defined as 


a 


2 






Z = 1 - - y2 _ 


(B.2) 


These variables are of course dynamical functions in general. Moreover, the characteristics 
of the coupled fluid, both at the background and perturbation levels, should be in general 
considered as functions of time, namely 


-|- [a (1 -b wq) — 1] Z 2 woa'^ -|- a [a (1 -b wq) — 1] Z 


(B.3) 


The derivative terms defined in (4.1), that are relevant for the perturbation evolution, for 
the coupling case (B.l) become 


X = —Kf3Z- 


a{l + Wo) 


(1 -b w) (0-2 -b Z) 

_ o I Za{l + WQ) 

~ (1 + w) (a^ + Z) ’ 


K^Z 

y = - 


Ct2 + Z’ 


z = k/3 


z 


Cj2 + Z 


r = — 


„ Za (1 -b Wq) [a (1 -b Wo) - 1] 

/-I . \ / O . r7\ 


(B.4) 


Hence, we can rewrite the evolution equation (4.11) completely in terms of the new variables 
(B.2), using (B.3),(B.4) above. This can be useful when elaborating the cosmological 
scenario numerically in terms of the background phase-space variables, but at this stage it 
might be not clear what we have gained. 

The benefit is obvious at the fixed points of the background. At these points the phase- 
space variables are constants, given by the four parameters of the model (a,/3, rco, A), and 
H scales as 


H = ——-—— , Weff = wo- {wo - 1) - (1 -b Wo) -b (1 -b Wo) {a-l)Z , 

3(1 + Weff) it - to) 

(B.5) 

with to an intergration constant. Therefore we can reduce the evolution equation (4.11) to 
the form 


5 + 


Cl 


+ + 


C 2 


t-to {t- to)' 


-S = k^c3 


t - to 


3(1+11 


1//1 


(B.6) 


where each of the three dimensionless coefficients q = Ci{a, (5, wo, A) is a constant given by 
the corresponding background fixed-point solution. Several different types of fixed points 
have been found for the model (B.l) and their stability properties have been analyzed in 
[43]. Given any of these solutions, one can immediately apply (B.6) and easily solve it to 
obtain the growth rate at a given scale as a function of the redshift. In general though, 
these expressions are rather complicated functions, and for simplicity we will not write 
them explicitly. 


^°Denoted a, x, y and z, respectively, in Ref. [43]. Here we use capital letters to avoid confusion with 
the definitions (4.1). We also used a slightly different notation from [43] in (B.l), where the parameters 
B and Vb are included in order to acquire the correct dimensions. 
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Let us only consider explicitly a special case for the sake of concreteness. Prom (B.3) 
we can confirm that the effective sound speed vanishes identically only if u^o = 0 and 
a = 1 (excluding the trivial case a = 0). This is of course what was expected as it is then 
that the Lagrangian contains only linear dependence upon the matter number density. 
Specializing to this case, we are left with the model parameters /3 and A. Actually, we have 
then reduced the model to a “conformally coupled” CDM, where p = m{^)n with now 
a non-trivial function m{4>) = 1 + {B (cf. (4.8),(4.16)). There is a well-known 
matter-scaling fixed point for the exponential quintessence for which X = Y = y^ 3/2/A and 
Z = 0. For this fixed point our equation (B.6) reduces to (setting to = 0 for convenience) 


? 4 j 2 / 3 , , 


(B.7) 


which is easily solved to yield the growth rate d log (5/(i log a = | — 72/A2 — 1^. In 

the limit A —)• oo this reduces to the usual matter dominated result d log 5/d log a —)• 1. 

To look at another hxed point for which the coupling is signihcant too, let us consider 


the exact background solution A = •\/1 , Y = 


_ ^ I Z+2p{l3-X) ^ ^ A(A-/3)-3 


(/3-A) 


- 1-1 


namely the 


conditionally stable fixed point F of [43]. At this fixed point, the equation (B.6) becomes 

(B.8) 


I ^ 4(4/3-A) j _ 2 [2/33A - 2/3^ (A^ - 12) + /3A - A^ + 3] ^ 


3Xt 


3\H‘^ 


It is again easy to solve this equation analytically. Now the resulting growth rate is however 
a more complicated function of the two parameters (the slopes of the coupling and the 
potential exponential functions). In order to present the results in a more transparent way, 
in Fig. 1 we plot the resulting growth-rate function as a function of the parameter /3, for 
four choices of the parameter A (note that physically the fixed point does not exist for all 
parameter values displayed in the hgure, namely it exists only for /3(/3 — A) > —3/2 [43]). 

B.2 Derivative couplings 

Here we adopt the model specified by 


/ = -^Ho/{Kn), 


(B.9) 


whose background phase space analysis was performed in [44] for the case with an ex¬ 
ponential potential. In analogy to our computations of the above subsection, where we 
used the previous background solutions for the algebraically coupled models, we can now 
employ the exact fixed-points solutions from [44] for the derivatively coupled models. Let 
us first look at the matter dominated fixed point p ^ (jP‘, V. At this fixed point Eq. (4.27) 
becomes 

S + —5 = -^ ^ „ J . (B.IO) 
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Figure 1. The growth rate for matter densities for the model (B.l) with a = 1, at the scalar-field 
dominated background solution (fixed point F of [43]). The x-axis is the coupling parameter [3 and 
we plot the growing solution for four choices of the potential parameter X as indicated in the legend. 


An analytic solution exists (with A and B dimensionless integration constants), namely: 
A 


5 = 


2F1 


Hot 
_ 1 

~ t 6 Js 


1 

1 

1 

CO 

2’ 

3’ 

6’ 

4 




+ Bmt^2Fi 


1 l.li 

3’2’ 6 ’ 4^ ° 




V2 

a + oscillations 


(B.ll) 


where 2 Fi{a, 6; c, z), is the Gaussian hypergeometric function and Jn{z) is the Bessel func¬ 
tion of the hrst kind of order n. In the second line of (B.ll) we have included only the 
approximate solution for the growing mode, while in the third line we have taken into 
account just the dominating small-argument limit of the Bessel function. As expected, we 
thus find again that to leading order dlog J/dlog a = 1, but there can occur slight modifica¬ 
tions as ^Hot begins to grow and become non-negligible. The details of the corresponding 
possible signatures in the matter power spectrum, and their impact to the constraints on 
the specific model, need to be computed by integrating the equations numerically. 

Finally, we will look at the same matter-scaling solution as for the previous model. At 
the hxed point B of Ref. [44], Eq. (4.27) becomes 




3t - eX^^Hot -t 4 (A2 - 3)" 


S = 


8 (A2 - 3)' 


3A4 -6X3(Hot + 4{X^-3f B 


(B.12) 

where, as explained in subsection 4.2.1, due to the choices ujq = = 0, / ~ 1/n, there is no 

scale-dependence. For completeness and transparency, we solve this equation numerically, 
and in Fig. 2 we present the resulting growth-rate evolution as a function of the scale factor, 
for six combinations of the model parameters A and As we observe, at high redshifts the 
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Figure 2. The growth rate for matter densities for the model (B.9) at the matter-scaling background 
solution (fixed point B of [44]) > a function of the scale factor, for six combinations of the model 
parameters X and f as indicated in the legend. The solutions begin from the constant growth rate 
given as the dominating solution to (B.7). 

solutions begin from the constant value given as the growing solution to Eq. (B.7). As the 
coupling becomes dynamically significant, the growth rate begins to evolve non-trivially, 
until the perturbations eventually freeze. However, the late evolution of the solutions should 
not be trusted as precisely when the coupling becomes significant, since the background 
solutions begin to deviate from the matter-scaling fixed point solution that was assumed 
in (B.12). Results from full evolution of the perturbations and systematic investigation of 
the parameter space of viable models, will hopefully be presented elsewhere. 

C Generalizations of the derivative couplings 

In this appendix we briefly consider the applicability of the scalar-fluid formalism within 
the context of so called disformal couplings [41, 42, 47, 59, 73]. In particular, we propose 
some possible generalizations of the scalar-fluid framework, in order to incorporate more 
general derivative-type interactions terms than the one in (2.21) in which we restricted to 
in this work. The disformal relation is a transformation of the metric of the form 

gpu = C{(]))g^u + . (C.l) 

The functions C and D could depend also upon the kinetic term of the scalar field, X = 
but for the purposes of the present discussion our simple points can be deduced 
as well from the less general starting point (C.l). 

Let us first consider just the scalar-fluid - formalism coupled in the disformal metric. 


- 24 - 





Then Eqs. (2.20),(2.21) generalize to 

p{h) + + J 

= X)p{n, s) + C((^)7-'(^, X)/(h, + J. (C.2) 

Here we have defined the function 7(0, X) as 7 = (l + ^Xj ^. This already suggests that 
the functions p and / in Eqs. (2.20),(2.21) could be considered to be functions of X as 
well^^. The number density in the disformal frame is given by 


n = 


C-l(<P)7(cl>,X) 



D{(l))gf,uJf^J^ 


(C.3) 


In practice it would be definitely easier to work in the disformal frame and obtain precisely 
the same sets of equations as in subsection 2.1, but with the appropriate tilde quantities, 
and then restore the metric g^u using (C.l). 

We can also follow the opposite approach, by starting from a disformal field theory 
and consider it in terms of the scalar-fluid formalism. As an example, we consider the 
Lagrangian of a test particle living in a disformal metric: 


C = moy/-gij,uXi^x^5^'^'> {x - x{X)) . 


(C.4) 


Varying the Lagrangian with respect to the metric g^u we obtain the stress energy tensor 
(2.24) of a pressureless perfect fluid, = pUfJJu, where we have made the identifications 


p = m(0) 

= 


I g^uxi^x 


g 


" ' 1 + ^ (<^,aC/“)' 


Ci^) 


5*-^^ {x — x(A)) 


X‘ 




y/-gapx°‘x>^ 


(C.5) 


We have denoted the time-dependent effective mass by m{4>) = and since 

the stress energy density is that of a dust-like component, we are allowed to make the 
straightforward interpretation that p = m(4>)n. We can then - at least formally - write 


£ = yf^p{(j),n) 


1 + 


D{4>) 

Ci4>) 




y/X^m{4>)n + 


D{^) m{<P) 


(C.6) 


The second form follows by naively using the relations of the fluid variables deduced in 
Section 2. Though this may not result in a completely equivalent theory in the scalar-fluid 
formalism, it clearly suggests the extension of the latter to allow for nonlinear derivative 
couplings between the fluid four-velocity and the gradient of the scalar field. One notes that 
the novel term appears naturally with the desired inverse proportionality to the number 
density (the metric determinant comes only to adjust the weight correctly). 

^^Note that the constraints are not affected since the metric decouples completely from J. However, one 
could formulate the theory in terms of instead, or consider = y/—gJ^, and thus the disformal relation 
(C.l) would result in quite non-trivial modifications. 
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